u=(t-tmin)/(tmax-tmin);
value=(1+sin(8*PI*p1*u+PI*p2))/2;

a=2*x+PI/2;
b=2*y+2*PI*p1;
c=2*z+2*PI*p2;
red=(1+sin(a))/2;
green=(1+sin(b))/2;
blue=(1+sin(c))/2;

x=20*(u-0.5);
y=20*(v-0.5);
r=(x*x+y*y)/8;
z=(1+p2)*6*cos(r+p1*PI)/(1+r);

red=(1+sin(2*PI*(u+p1)))/2;
green=v;

dx=-x;
dy=-y;
dz=-z;

red=(1+x)/2;
green=(1+y)/2;
blue=(1+z)/2;

nu=u*(abs(p1)*29+1);
nv=v*(abs(p2)*29+1);
a=nu-floor(nu);
b=nv-floor(nv);
red=(a>0.5);
green=(a>0.5)*(b>0.5);
blue=(b>0.5);